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I.  SUMMARY 


Yhis  is  a  study  of  vibration  control  for  large  space  structures.  Advantage  is 
taken  of  the  limiting  performance  characteristics  of  dynamic  systems.  This  approach 
permits  large  problems  with  constraints  to  be  analyzed.  A  modal  formulation  for  the 
limiting  performance  was  developed  in  order  to  enhance  the  applicability  of  limiting 
performance  to  large  structural  systems.  One  effort  to  develop  an  optimal  control 
system  is  based  on  the  limiting  performance  approach  in  combination  with 
classical/optimal  control  theory.  A  limiting-performance/minimum-time  solution  was 
formulated  to  achieve  the  goal  of  rapid  suppression  of  disturbances.  Classical/optimal 
control  studies  show  that  a  position  loop  might  be  useful  in  taking  care  of  constraint 
controllers,  such  as  proof— mass  dampers.  Finally,  to  derive  feedback  control  law 
based  on  the  limiting  performance  characteristics,  parameter  identification  technique  has 
been  under  investigation. 


II.  RESEARCH  OBJECTIVES 


1 .  Modal  Limiting  Performance  Formulation 

Instead  of  evaluating  Duhamel  integrals  to  implement  a  linear 
programming  formulation  of  the  limiting  performance  problem,  a  modal 
approach  was  to  be  studied.  The  focus  of  the  study  was  to  develop  a 
procedure  for  systemizing  the  coupling  of  vibration  controllers  to  an 
existing  (modal)  structural  model.  A  study  of  vibration  control  of 
structural  systems  was  to  be  performed. 

2.  Limiting  Performance  Control  of  Large  Space  Structures 

The  limiting  performance  formulation  was  to  be  enhanced  and 
expanded  to  treat  control  problems  of  large  space  structures  (LSS)  with 
proof-mass  actuators.  To  take  the  characteristics  of  the  LSS  problem 
into  consideration,  a  minimum-time  solution  ,  which  was  intended  to 
dissipate  the  undesired  vibration  rapidly,  was  considered.  Methods  to 
incorporate  a  minimum  time  solution  into  the  limiting  performance 
formulation  were  of  major  concern. 

3.  Inertia  Dampers 

In  addition  to  the  limiting  performance  approach,  inertia  damper 
design  problem  was  to  be  treated  using  the  classical  or  optimal  control 
laws.  It  was  expected  that  the  success  of  the  control  laws  could  be 
measured  by  the  limiting  performance  control  law.  If  the  resulting 
control  laws  give  characteristics  which  are  close  to  those  of  the 
limiting  performance,  the  control  systems  may  be  considered  to  be 
optimal.  Since  the  inertia  dampers  were  subject  to  constraints,  a  major 
consideration  was  related  to  the  methodology  for  incorporating  these 
constraints  into  the  classical  or  optimal  control  laws. 
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•4 .  Identification  of  Conrro]  laws 


In  a  preliminary  study,  identification  of  control  laws  based  on  the 
limiting  performance  characteristics  was  to  be  studied.  This  technique 
can  handle  the  problem  efficiently  because  once  the  limiting  preformance 
characteristics  are  known,  the  identification  process  can  be  separated 
from  the  system  dynamics. 


III.  RESEARCH  PROGRESS  AND  STATUS 


I .  Background  Information 

Formulation  oF  the  Linear  Programming  Problem 

A  linear  vibrating  system  with  n  degrees  of  freedom  subjected  to 
arbitrary  excitation  f(t)  and  control  force  u(t)  is  described  by  the 
equat ions 


Mx  +  Cx  +  Kx  +  Vu  =  f(t) 


(1) 


where  x  is  the  displacement  vector,  M,  C,  K  are  the  n  x  n  mass,  damping, 
and  stiffness  matrices,  respectively;  the  n  x  j  coefficient  matrix  V, 

associated  with  the  j  control  forces  u^  .  .  u^.,  places  these 

forces  in  the  apppropriate  rows  in  Eq .  (1).  In  general,  the  measure  of 
performance  may  be  a  linear  combination  of  displacements,  velocities, 
accelerations,  control  force  components,  and  external  forces.  To  define 
such  a  performance  index,  let 


h  = 


PjX  + 


P2x  ♦ 


P3x  ♦ 


p4y 


+  pr  f 


(2) 


where  the  P^  are  prescribed  coefficient  matrices.  It  is  clear  that  h 
may  be  considered  to  be  explicitly  dependent  on  time  t  and  the  control 
forces  u. 


h  =  h  ( t  ,  u  ) 


(3) 
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The  performance  index  +  is  then  written  in  terms  of  h  as 

=  max  max  |  h  ( t .  u )  |  (4) 

i  t  i 

where  i  varies  over  the  rows  of  h  and  t  over  the  time  interval  of 
concern.  With  these  definitions,  the  optimization  problem  becomes 
minimization  of  >//  with  respect  to  the  control  forces  u^,  .  u  ^ . 

For  any  system  governed  by  Eq.  (1).  the  system  responses  are  linear 
functions  of  u.  It  is  convenient  at  this  point  to  discretize  u  in  the 
time  domain.  Thus,  if  the  time  interval  of  interest  is  0  <  t  <  Tj.  and 
this  interval  is  represented  by  N  discrete  instances  of  time,  then  u 
takes  on  N-l  discrete  values.  This  piecewise  constant  discretized 
function  is  denoted  u.  Then,  the  discretized  form  of  Eq .  (2)  becomes 


h=Wu+g(t)  ( 5 ) 

where  g  is  an  explicit,  known  function  of  time  and  W  is  obtained  by 
solving  Eq.  (1)  for  the  response  variables  in  terms  of  u,  discretizing 
the  results,  and  replacing  the  first  four  terms  of  Eq .  (2)  with  these 
resul ts . 

To  put  the  optimization  procedure  in  standard  linear  programming 
form,  define 


b 


u(l) 

u(2) 


(C) 


L  u(N-l)  J 


where  u(k)  means  the  vector  u(t)  evaluated  at  the  k-th  subinterval  of 
time.  Thus,  if  the  time  discretization  is  uniform  with  each  subinterval 
h  seconds  long,  u(k)  =  u((k-l)h).  Next,  let 


z  = 


and 


c  =  [1  0  ...  3] 


(7) 


(S) 


Tlien  the  linear  programming  problem  is  to  minimize 


4*  =  c  z 


subject  to  the  constraints 


Hz  <  b 


The  constraints  (10)  always  include  the  conditions 


(9) 


(10) 


0 


A 


A 


'P  >  |h.  (t  ,u)  | 


(11) 


Thus,  minimizing  over  the  u  such  that  conditions  (11)  hold  makes  ■ p  the 
least  upper  bound  of  the  set  {|h.(t.  u)|}  over  time  t  and  control  forces 
u.  In  addition  to  the  constraints  (11).  Eq.  (10)  may  also  contain 
bounding  constraints  on  response  variables  or  control  force  magnitudes. 

2 .  Modal  Formulation  F 1 1 

A  modal  formulation  of  the  limiting  performance  problem  has  been 
proposed  as  one  of  the  tasks.  Instead  of  evaluating  Duhamel  intergrals 
to  accomplish  the  linear  programming  formulation,  an  alternative 
approach,  a  modal  approach  was  explored.  This  approach  is 
computationally  convenient  because  it  simplifies  the  steps  necessary 
before  linear  programming  optimization  can  be  initiated.  Also,  modal 
trunctation  can  be  incorporated  easily  to  handle  large  finite  element 
models  such  as  large  space  structure  control  problems.  Furthermore,  the 
practicality  of  determining  limiting  performance  characteristics  is 
enhanced  because  modal  properties  may  be  experimentally  available.  This 
modal  formulation  was  applied  to  find  the  limiting  performance 
characteristic  of  large  structures,  perhaps  formed  of  substructures 
subject  to  transient  distrubances.  Modal  formulations  were  dclevoped 
for  both  systems  with  imbedded  control  forces  and  systems  with  generic 
control  force  connections  between  substructures. 


Undamped  Systems 


for  an  undamped  system. 


mode  shapes  <£ 

k 


and  mode  frequencies  may 

K 


be  found  from  the  eigenvalue  problem 


(-wjW)  =  o  (12) 

If  mode  shapes  are  dimensional ized  such  that 

$TMf  =  I  (13) 

$TK<I>  =  n2  =  diag(w2  Ji . u2)  (H) 

-  -  ~  '‘12  ny  v 


where 


?  =  [  ±2  -n  ]  (15) 

then  tlie  solution  to  Eq.  (1)  is  given  by 

n  n 

x(t)  =  k2j  cosc^t  4k«kM  x(0)  +k“15lnU’kt\»1.M  x(0) 

uk 

n  t 

+  kll  -k  -  ^  -k  U(T)_  VH  ( T )  ]  sinc.'k(  t-r)dr  (16) 

w,  0 

K 

No  conditions  have  been  placed  on  the  control  forces  which  replace  the 
design  elements:  they  may  be  linear  or  nonlinear.  In  the  case  where 
control  forces  have  replaced  passive  elements  such  as  springs  or 
dashpots,  it  should  be  clear  that  the  mode'  shapes  and  frequencies  in  Iq. 
(16)  are  those  for  a  system  without  these  elements  included.  The 
contributions  to  the  displacement  coordinate  solution  (16)  from  the 


S 


replaced  elements  are  contained  in  the  control  forces  which  have 
replaced  them. 

The  advantage  of  using  a  modal  formulation  for  the  response 
variables  is  especially  apparent  when  dealing  with  systems  with  control 
force  connections  between  components.  Mode  frequencies  in  Eq .  (16) 
are  identically  those  of  the  components  and  mode  shapes  ^  are  readily 
available  from 

[  *1  -2  .  .  . -n  ] 

(17) 

where  is  a  matrix  whose  vectors  are  the  modeshapes  for  the  i-th 
component.  Modal  properties  for  each  component  may  be  obtained 
experimentally  or  determined  computationally  with  general  purpose  (FE) 
computer  programs. 

If  control  forces  are  discretized  in  some  fashion  then  the 
discretized  response  at  time  t=mT  is  obtained  from  Eq .  (16) 

x ( mT )  =  r ( mT )  (IS) 

n  .  T  m- 1 

-  2  -2  2  u  (  i  T )  [cos  (m-i-l)w.T  -  cos(m-i  )w  ' T] 

k,  (x),  K  K  .  ~  K  K 

=  1  k  i=0 

and  r(mT)  represents  contributions  from  the  disturbance  forces  f  and 
initial  conditions.  Equation  (16)  may  be  differentiated  and  discretized 
to  obtain  velocity  and  acceleration  expressions 


<) 


(ID) 


x  (  mT  )  =  r  ( mT ) 


n  .  in- 1 

-  2  -  <#>.  <#> 1  V  2  u(iT)  [sin(m-i)u,T  -  sin(m-i-l )w,  T] 

.  .  (J.  — K~K  -  n  ~  ^  K 

k=l  k  1=0 


x(mT)  =  r(mT)  -  2  <£ ,0 ,Vu(mT) 

k=l  R  k  ~ 


n  y  ml 

-  2  2  u(iT)  [cos(m-i  )w.  T  -  cos(m-i-l  )u,T] 

k=rk_k  i=o  k  k 


In  matrix  notation,  Eqs.  (IS)  to  (20)  have  the  form 


x(mT)  =  r(mT)  +  Qu 
x(mT)  =  r(mT)  +  Ru 
x(mT)  =  r(mT)  +  Su 


where  Q,  R,  and  S  are  dependent  upon  modal  properties,  V,  and  the 
discrete  time  interval  chosen.  With  the  use  of  these  relations,  the 
vector  h  as  a  function  of  control  forces  and  time  is  given  by 


h  =  (P^  +  P2R  +  P3Q  +  P4  )  u  +  P5f 
+  P j  r ( mT )  +  P2r(mT)  +  P3r(mT) 


Once  response  variables  and  h  are  expressed  as  functions  of  u,  the 
performance  index  and  all  constraints  may  be  written  in  terms  of  the 


vector  7.  for  each  discrete  instance  in  time. 
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A  similar  analysis  may  be  obtained  for  damped  systems.  When 


viscous  damping  is  present,  mode  shapes  4  and  eigenvalues  X^  may  be 
found  from  the  eigenvalue  problem 

(Xk  M  +  AkC  +  =  -  1  <  k  <  2n  (25) 

If  mode  shapes  are  dimensional i zed  such  that 

i  (2AkM  +  C)4  =  1  (26) 

then  the  solution  to  Eq.  (1)  is  given  by 

2s  X,  t  t  -X,  t  _ 

x(t)  =  r(t)  -  2  4  e  /  e  ^  Vu(T)dT 

k=l  “K  0  k 

n+S  ,  C  Mt_T)  T 

-2  2/  Re(e  *>.'iOVu(r)dT  (27) 

k=2s+l  0  “k  K 

where  it  is  assumed  that  the  first  2s  eigenvalues  are  real  and  X 

&  n+s+1 

...  ,  X_  are  the  complex  conjugates  of  X0  . .  X  .  r(t)  denotes 

2s+l  n+s 

the  part  of  the  response  independent  of  the  control  forces. 

If  control  forces  are  discretized  in  the  same  fashion  as  for  the 
undamped  case,  the  discretized  response  at  time  t=mT  is  obtained  from 


T 


2s  ^k*"  ~  m-1  -A,(i+1)T  -X,  iT 

x ( m T )  =  r(mT)  +  2  e  2  u(iT)  [e  -  e  ] 

k=l  X,  i=0 

k 

n+s  ,  T  m-1  X,(m-(i+l))T  X,(m-i)T 

+  22  Re{i  +  1  V  2  u(iT)  [e  -  e  ]}  (28) 

k=2s+l  k  i=0 

Velocity  and  acceleration  expressions  may  be  found  by  differentiating 
Eq.  (27)  and  discretizing  the  result.  The  remaining  steps  necessary  to 
define  the  linear  programming  problem  are  not  presented  since  they  are 
identical  to  those  followed  for  the  undamped  case. 

Appl ication 

A  three  DOF  system  was  treated  to  demonstrate  the  concepts 
presented  in  modal  formulation. 

3 .  Limiting-Performance/Minimum-Time  Formulation  f 31 

A  minimum  time  solution  has  been  superimposed  on  the  conventional 
limiting  performance  response  to  achieve  the  rapid  suppression  of  the 
disturbances  in  the  minimum  time.  Since  the  min-max  norm  of  the 
limiting  performance  gives  a  unique  solution  only  until  the  peak  value 
of  the  performance  index  is  achieved,  an  additional  measure  of 
performance  is  desired  to  obtain  a  unique  solution  after  the  peak  value. 
Two  different  approaches  were  studied  to  accomplish  the 
1 imi t ing-per f ormance/minimum-t ime  solution.  One  approach  uses 
additional  constraints  and  the  other  is  based  on  the  performance  index. 
The  two  methods  were  applied  to  the  control  of  a  simple  model  of  a 
conti lever  beam  with  a  proof-mass  damper  and  the  two  methods  led  to  the 
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T 


identical  minimum  time  solutions. 

Problem  Statement 

A  linear  vibrating  system  with  n  degrees  of  freedom  subject  to 
arbitrary  external  excitations  £(t)  and  control  forces  u(t)  is  expressed 
in  the  first  order  system  of  differential  equations 

s(t)  =  As( t )  +  Bu( t )  +  Cf(t)  (2D) 

where  s(t)  is  an  n-dimensional  state  vector.  A,  B,  and  C  are  n  x  n,  n  x 
nu  and  n  x  nf  constant  coefficient  matrices.  The  quantities  nu  and  nf 
are  the  number  of  control  forces  and  excitations,  respectively. 

Constraints  are  imposed  on  the  dynamic  system  under  study.  The 
format  of  the  constraints  is 

<  QjS  +  QgU  +  Q3I  <  Yy  for  t  <  t  <  tf  (30) 

where  y^  and  y^  are  nc-dimensional  lower  and  upper  constraint  vectors, 

Qj ,  Qg  and  are  nc  x  n,  nc  x  nu  and  nc  x  nf  constant  coefficient 
matrices,  and  t  and  t^.  are  given  initial  and  final  time. 

The  problem  is  to  find  an  optimal  control  u(t)  which  will  transfer 

an  initial  state  s(t  )  =  s  to  a  desired  final  s'tate  s(tr)  =  sr  in  the 
-v  o'  _o  ~v  f'  — f 

minimum  time  while  extremizing  a  given  performance  index  of  the  form 

maX  rp  rp  rp 

Minimize  J  =  { tQ  <  t  <  t^.  jp^s  +  ££2  +  P3- 1 }  (31) 
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where  £2  an<^  £3  are  given  n.  nu  and  nf  constant  coefficient  vectors. 

Formulation  Using  Constraints 

In  addition  to  the  constraints  in  Eq.  (30),  an  additional 
constraint  set,  referred  to  as  steady-state  constraints,  is  imposed  on 
the  time  from  the  assumed  minimum  time  { t  ^ )  to  the  final  time  (t^.). 

%L  -  Q1S  -  +  Q2S  y  +  Q3S  -  -  ySU  for  Ct  -  C  -  lf 


where  Yg^  and  Yg^  are  ncs-dimensional  coefficient  vectors  representing 
lower  and  upper  bounds  of  steady-state  constraints;  Qjg.  Q^S’  Qgg 

are  ncs  x  n,  ncs  x  nu,  and  ncs  x  nf  coefficient  matrices.  The 
steady-state  constraints  represent  the  desired  final  states  of  the 
system. 


To  place  the  optimization  procedure  into  the  standard  linear 
programming  form,  the  system  in  Eq .  (29)  is  discretized  using  uniform 
time  intervals  to  obtain  a  set  of  state  difference  equations 

s(k+l )  =  Gs(k)  +  H[Bu(k)  +  Cf(k)]  (33) 


where 


s(k)  =  state  vector  at  time  t  =  t^ 

u(k),  £(k)  =  control  and  external  excitation  vector  at  t  =  t,  , 

K 

assumed  to  be  constant  over  the  interval  t,  <  t  <  t,  , 

k  ~  k+1 


H 


H  =  /h  eA(h'T)dT 

o 

h  =  time  interval  =  t^+^ 
k  =  1,  2 . N-l 


t 


k 


The  state  vector,  at  any  time  t  =  t^,  can  be  expressed  as  a 

function  of  the  initial  state  s(l),  the  control  history  u( 1 ) .  u(2) . 

u(N-l)  and  the  external  excitation  f(l).  £(2) . £(N-1). 

,  k- 1 

s(k+l)  =  GK  s(l)  +  2  Gk  J  H[Bu(j)  +  Cf(j)] 

j  =  l 

(3-1) 

+  H[Bu(k)  +  Cf(k)] 

k=  1,  2 . N-l 


The  constraints  in  Eq.  (30)  and  (32)  are  discretized  similarly 

YL(k)  <  Q1s(k)  +  Q2ij(k)  +  Q3£(k)  <  y  (k)  for  k  =  1.2 . N-l  (35) 

YSL(k)  1  <>1S2(k>  +  Q2Sy(k)  +Q3S-(k)  1  YSU(k)  (36) 

for  k  =  NT.NT+1 . N-l 


where  NT  is  discretized  assumed  minimum  time. 

The  objective  function  of  Eq.  (31),  which  reflects  the  min-max 
norm,  is  discretized  and  converted  into  a  constraint  set.  Since  J  is 

X  T  T 

the  maximum  value  of  |p^s  +  p^u  +  I 
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Ie[s(k)  +  ( k )  +  E3I ( k )  I  1  -I 

-J  +  CBjs(k)  +  e^u(k)  +  EgKk)]  <  0 
-j  -  [E|s(k)  +  +  E^K^)]  <  0 


Def ine 


(V/) 


(3S) 


(39) 


subject  to  the  constraints 


Hz  <  b  (43) 


where  II  and  b  represent  constraints  of  Fqs.  (3a),  (3C>).  and  (3S). 


1(> 


VtTi  i  1  o  extremi/ing  the  given  performance  index,  the  linear 

programming  routine  will  cheek  the  feasibility  of  the  solution  with  the 

given  initial  states,  external  forces,  and  the  two  sets  of  constraints. 

The  minimum  time  is  the  smallest  time  which  will  make  the  solution 

feasible  with  the  given  conditions.  If  the  response  beltaves 

monotonical  ly ,  i.e.,  increases  or  decreases  without  fluctuation,  t^.  =  t^ 

can  be  used.  However,  if  the  response  varies  sinusoidally  a  very 

careful  selection  of  t^.  is  required.  For  the  sinusoidal  response,  t  j.  - 

t  must  be  at  least  a  half  period  of  the  response.  The  search  for  the 

smallest  feasible  time,  t  ,  can  be  done  efficiently  using  the 

mm 

bisection  method. 


Formulation  Using  Performance  Index 

To  achieve  alternative  1 imi t ing-porformance/minimum-t ime  solution, 
the  performance  index  is  given  in  F.q .  (31)  is  modified.  Two  sets  of 
performance  indices  are  considered.  One  set  of  them,  referred  to  as  the 
transient  performance  index,  is  given  by 


Jt  "  r  < 
o  - 


max 


ax  I  T  T  i  .  | 

t  <  t  'El5  +  02>J  +  03f  I 


(11) 


where  t  is  the  time  limit  for  the  transient  period.  The  other  set, 
referred  to  as  the  steady-state  performance  index  is  defined  as 


.  max  I  T  T  Tr  | 

Js=  t  <  t  <  tf  ^  +  +  e.*£l 

t  -  f 


(T-) 


Now,  the  "global"  performance  index  is  defined  by 
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and 


J  =  Jt  +  Js 


(■Hi) 


Note  that  the  vectors  p^ •  and  O3  are  not  changed  in  Eqs.  (  )  and 

(•15). 


Consider  now  the  linear  programming  formulation.  The  objective 
functions  of  Eqs.  (-H)  and  (15)  which  reflect  the  min-max  norm  can  be 
converted  into  a  constant  set  as  follows: 


Ipf®  +  P2-  +  P3I I  1  Jt 

and  (-17) 

bjs  +  4-1  +  D3II  <  Js 

To  place  this  optimization  problem  into  a  standard  linear  programming 
form ,  def ine 


z  '  = 


(•IS) 


where  u  =  [u( 1  )  u(2)  ...  u(N  -  1)]^ 


(-10) 


and 

c'T  =  [1  1  0  ....  0] 


Then  the  linear  programming  problem  is  to  minimize 


(50) 


(51) 


IS 


subject  to  the  constraints 


H'z '  <  b'  (52) 

where  H'and  b‘ represent  constraints  of  Eos.  (30)  and  (17). 

The  minimum  time  (t  )  is  the  smallest  time  which  will  make  the 
mm 

global  performance  index  of  Eq .  (T6)  stay  within  a  desired  value.  Since 

the  performance  index  can  be  computed  for  each  iteration,  an 

interpolation  method  such  as  the  secant  method  can  be  employed  to  find 

t  ef  f icient iy . 
min  J 

T .  Suhoptjmal  Feedback  Control  for  the  I’roof-Mass  Actuator  Using 
a  Eos i  t  ion-l  oop  T2,  3) 

In  addition  to  the  various  limiting  performance  approaches,  the 
problem  of  the  optimal  design  of  the  proof-mass  dampers  in  a  large 
structural  system  using  the  linear  control  laws  and  the  optimal  control 
theory  was  studied.  Initially,  linear  control  laws  to  control  the 
proof-mass  dampers  were  investigated.  The  results  showed  that,  due  to 
the  constraints  of  the  proof-mass  dampers,  although  adequate  damping 
could  be  achieved  at  high  frequencies,  very  little  damping  could  be 
obtained  at  frequencies  of  one  Hz  or  less.  While  maintaining  the 
adequate  damping  at  higher  frequencies,  to  improve  damping  at  lower 
frequencies,  both  the  limiting  performance  technique  and  optimal  control 
technique  were  i nves t iga t ed  and  compared.  1  he  limiting- 
performance/mi  nimum- t ime  solutions  were  obtained  for  the  control  of  a 


11) 


1 


I 

I 

proof-mass  damper  attached  to  the  end  of  a  c  kimped- f ree  beam  which 
represented  a  small  scale  preliminary  model  of  large  1 lexible  space 
structures  and  used  as  a  measure  of  the  success  of  the  optimal  control 
law.  The  system  was  composed  of  a  proof-mass  actuator,  a 
servo-accelerometer  attached  to  the  free  end  of  the  beam  and  a 
proximeter  to  measure  the  position  of  the  proof-mass  along  its  track. 

Using  the  optimal  control  theory  and  the  full  state  feedback,  a  control 
law  which  was  simple  but  stable  both  in  the  linear  and  nonlinear  region 
was  obtained.  To  take  the  constraints  into  consideration,  a  position 
loop  was  added  to  the  system  and  the  constraints  of  rattlespace  and 
control  force  of  a  proof-mass  actuator  were  able  to  be  treated  by  only 
limiting  the  desired  relative  position  of  the  proof-mass.  The  resulting 
suboptimal  control  law  showed  adequately  damped  response  at  low 
frequencies  which  was  very  close  to  tliat  of  limiting-performance/ 
minimum-time  control. 

A  Mathematical  Model 

A  perfect  model  is  rarely  obtained  for  any  physical  process.  One 
of  the  most  prominent  sources  of  modeling  error  in  large  flexible  space 
structures  is  the  deletion  of  modes  in  the  formation  of  the  design 
model.  The  low-frequency  modes,  which  are  more  accurately  known,  are 
retained  and  the  high-frequency  modes  are  deleted.  Assume  a  perfect 
actuator  and  consider  only  the  first  mode  of  the  beam  to  be  perfectly 
known,  identified,  and  observed.  The  residual  unknown  modes  are 
represented  by  the  second  mode.  This  mode  usually  has  eigenf requency 
that  is  not  far  from  the  closed-loop  bandwidth  and  must  be  considered  in 
the  control  system  design.  Modern  control  theory  is  used  to  obtain 
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optimal  first  mode  performance  and  a  stable  second  mode.  It  is  assumed 
that  higher  residual  modes  have  e igenf ' recjunc i es  above  the  main  control 
system  bandwidth.  They  can  be  damped  by  passive  dampers  and  do  not 
interfere  with  the  low-frequency  control  law.  Assume  that  the  plant's 
variations  in  time  tire  much  slower  than  its  characteristic  time 
constant,  so  it  can  he  treated  as  a  time  invariant  system.  A  simplified 
model  for  the  first  two  modes  of  the  beam's  free  end  is  treated.  The 
system  is  subjected  to  the  following  control  force  (u)  and  rat  t lespace 
(d)  constraints 

-u  <.  u  <  u 

max  -  -  max  (5.3) 

~d  <  d  <  d 

max  -  -  max  (51) 

h1  :hn'i  t  i  ,r*u  1  Con  ;  ml  ( 'ops  jder  i f  bus  t  ra  i  n  t  s 

It  is  apparent  from  the  results  of  the  minimum  time  solution  that 
the  prool -mass  relative'  motion  alon.;  its  finite  length  track  is  rather 
simple.  lhis  kitnl  of  motion  can  be'  implemented  by  usiryj;  a  direct 
position  control.  I  he  system  is  expanded  by  considering  a  position  loop 
added  t e>  the  plant.  The?  input  to  the  system  is  the  desired  relative 
position  of  the  proof-mass  d  ,  where  d  is  its  actual  relative  position. 
Int roduct it>n  of  a  position  loop  into  the  plant  and  the  proper  choice  of 
its  trains  make  it  possible  to  reduce  the  problem  to  a  constrained  input 
proh  1  i'tn . 


for  the  performance  index 


the  opt  i  1 1  control  can  he  treated  through  the  calculus  of  variations 
and  Hamilton's  principle.  The  Hamiltonian  is  defined  by 


where  ,\  is  the  costate  vector. 

In  the  case  of  constrained  input,  it  was  shown  by  Pontryngin  et  nl. 

that  "the  Hamiltonian  must  be  minimized  over  all  admissible  control  for 

optimal  values  of  the  state  and  costate".  In  the  case  under 

consideration,  the  Hamiltonian  is  a  quadratic  function  of  the  control 

d  There  will  he  only  one  location  which  satisfies 
c 


if  i! 


c: 


(57) 


This  point  will  be  denoted  as  cl  ,d  derived  from  the  solution  of  the 

c 

Riccati  equation  .  For  constrained  input.  Pont ryagin ' s  Minimum 
Principle  leads  to  the  following  solution 


d 

max 

for 

3  > 

c 

d 

FTHIX 

d  - 

c 

cl 

c 

for 

|3  1  < 

c 

d 

nvix 

(58) 

-d 

nv  ix 

for 

3  < 

c 

-d 

rrt  ix 

Interest  is  in  the'  steady  state  behavior  of  the  process  where  the 


optimal  feedback  gain  matrix  is  obtained  from  the  solution  of  the 
Algebraic  Hiccaii  Kquation  to  minimize  the  following  quadratic 
performance  index 


J  -  0.5 


[ Q  s  +  R“  d  ]  dt 
-  c 


(50) 


The  resulting  linear  position  command  is  represented  by 


a 


-  ( K  xt+K  x ,  +  K  ,  d  +  K  d) 
v  x  v  d  w  ’ 


(60) 


5 .  System  Identification  of  Snhoptimnl  Feedback  Control  Parameters 
Rased  on  the  limiting  Performance  Characteristics 


Since  the  limiting  performance  gives  the  best  possible  or 
"limiting"  response  of  a  system,  it  would  appear  to  be  reasonable  to 
base  a  control  system  on  the  limiting  performance  characteristics. 
However,  due  to  uncertainties  in  control  problems,  open  loop  control 
such  as  the  limiting  performance  control  may  not  be  applicable  in 
practice,  unless  real-time  computiu.,  power  for  the  limiting  performance 
is  available.  To  overcome  this  difficulty,  parameter  identification  to 
find  suboptimal  feedback  control  laws  based  on  the  limiting  performance 
characteristics  is  under  study. 

Consider  a  linear,  time  invariant  dynamic  system  represented  by  a 
set  of  difference  equations 

s(k  +  1)  =  Gs(k)  +  HBu(k)  (61) 


JL 
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From  the  limiting  performance  characteristics,  the  optimal  time 
responses  s  (k)  and  u  (k)  are  obtained. 

Consider  a  linear  controller. 

u(k)  =  Ks(k)  (62) 

where  K  is  an  m  x  n  feedback  gain  matrix. 

Since  m  controllers  are  considered  and  optimal  control  forces  are 
available  for  each  controller,  it  is  possible  to  proceed  controller  by 
controller.  Also,  the  system  is  assumed  to  be  subject  to  the 
constraints.  For  controller  i,  a  suboptimal  linear  control  law  to  be 
determined  is  described  by 

Umax.  for  Us.(k)  >  Umax , 

i  1  l 

l's.(k)  =  <  Us*(k)  for  |Us*(k)|  <  Unutx.  (63) 

for  Us^(k)  <  Umin. 

where 

USj(k)  =  suboptimal  control  force  for  controller  i 

U-Sj(k)  -  k.s*(k)  (61) 

.  th 


Umin . 

l 


k. 
- 1 


row  of  K  matrix 


Once  the  optimal  limiting  performance  characteristics  are  known,  it 
is  not  necessary  to  perform  a  complete  structural  analysis  for  this 
identification  procedure.  Thus,  this  identification  method  handles 
problems  efficiently. 

It  is  expected  that  this  pramether  identification  technique  will 
give  a  stable  feedback  control  law  for  a  large  structural  system.  To 
ensure  stability  the  Routh-Hurwi tz  criterion  is  considered  to  be 
incorporated  as  constraints  when  a  least  squares  fit  is  performed.  An 
efficient  computational  scheme  to  achieve  a  stable  feedback  control 
based  on  the  limiting  performance  characteristics  is  under  study. 
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